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ABSTRACT 

The long gamma-ray burst GRB 060714 was observed to exhibit a series of 
five X-ray flares beginning ~ 70 s after the burst trigger To and continuing until 
~ T + 200 s. The first two flares were detected by the Burst Alert Telescope 
(BAT) on the Swift satellite, before Swift had slewed to the burst location, while 
the last three flares were strongly detected by the X-Ray Telescope (XRT) but 
only weakly detected by the BAT. This burst provides an unusual opportunity 
to track a complete sequence of flares over a wide energy range. The flares were 
very similar in their light curve morphology, showing power-law rise and fall 
components, and in most cases significant sub- structure. The flares also showed 
strong evolution with time, both spectrally and temporally. The small time scale 
and large amplitude variability observed are incompatible with an external shock 
origin for the flares, and support instead late time sporadic activity either of the 
central source or of localized dissipation events within the outflow. We show 
that the flares in GRB 060714 cannot be the result of internal shocks in which 
the contrast in the Lorentz factor of the colliding shells is very small, and that 
this mechanism faces serious difficulties in most Swift GRBs. The morphological 
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similarity of the flares and the prompt emission and the gradual and continual 
evolution of the flares with time makes it difficult and arbitrary to draw a dividing 
line between the prompt emission and the flares. 

Subject headings: gamma rays: bursts 



Introduction 



One of the most surprising findings of the studies of gam ma-ray bursts (GRBs) made with the 
Swift Gamma- Ray Burst Explorer (IGehrels et al.l 120041 ) is that nearly half of all bursts show 
flares, or la rge short-lived increas e s in emission at times after the initial pr ompt emission has 
died away (IBurrows et al.l l2005bl ; iNousek et al.l 120061 ; lO'Brien et al.l 120061 ) . These flares are 
superimposed on either the rapidly decaying tail of the prompt emission, or the very slowly 
decaying phase of the early afterglow, and can involve flux increases of as much as three 
orders of magnitude. Some bursts have single flares, although most bursts with flares have 
multiple flaring episodes. Most flares are at early times, t < 10 3 s, although strong flares 
can occur as late as > 10 4 s after the onset of the burst. 

The origin of GRB flares is still an open question. Evidence is mounting, however, that 
the origin of these flares is similar to that of the prompt GRB emission (i.e. either internal 
shocks or some other well localized dissipation process within the ultra-relativistic outflow), 
rather than to that of the afterglow emission (i.e. the external shock going into the ambient 
medium) . This evidence includes the multiplicity of flares; their sharp time structure: rapid 
rise and decay and sub-peak structure within the flares; the large increase in flux during the 
flare; and the hard to soft spectral evolution of the flares, which is similar to the spectral 
evolution found in the prompt emission. 



In the context of the internal shocks model (IRees fc Meszarosl Il994l ). flares may be 
caused by late-time collisions of shells of relativistic material that are produced by the 
central engine with varying Lorentz factors. These can occur either by late time sporadic 
activity of the central source, or by a s mall relative velocity between shells that were ejected 
during the prompt GRB emission (e.g.. IBurrows et al.l 120071 ) . 



Very few GRBs have a sequenc e of multiple flares bright enough to be studied in detail. 
GRB 060714 (IKrimm et al.l |2006af) showed a series of five flares starting at ~ 70 s after 



the start of the burst. The first t 



iree f 



Telescope (BAT; iBarthelmv et al 



ares were clearly detected by the Swift Burst Alert 
2005 aj) and the last three were se en as strong flares by the 



Swift X-Ray Telescope (XRT; iHill et al.ll2004l : IBurrows et al.ll2005al ) and as weak flares by the 
BAT. The first two flares occurred while the spacecraft was slewing to or settling at the burst 
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location, so they were not observable by the XRT. This burst provides a rare opportunity to 
study a rapid sequence of flares across a large energy range. The five flares show evidence of 
hard to soft spectral evolution as the flares progress and also strong similarities, in particular 
sharply resolved temporal features and large flux increases in each flare. We show that these 
results are inconsistent with an external shock (i.e. afterglow) origin for the flares, and 
suggest instead a late-time and lower energy continuation of the prompt emission, either 
due to late time intermittent activity of the central source or via well localized spasmodic 
dissipation events within the outflow. 

In fact these results suggest that it may no longer be possible to draw a clear distinction 
between what have traditionally been called the prompt emission and the X-ray flares. His- 
torically the prompt emission was that detected above ~ 20 keV, and generally considered 
to be due to activity of the central engine. X-ray flares are usually detected most strongly 
at lower energies, so they were considered a completely separate phenomenon. Since we 
show here that (a) the flares of GRB 060714 are very likely of common origin to the earliest 
emission from the burst, (b) the flares are detected above 20 keV, and (c) the flares show a 
gradual and continual evolution linking them to the prompt emission, it is quite reasonable 
to consider the flares a lower energy continuation of the same phenomenon as the prompt 
emission. However to be consistent with earlier work we do use the term "prompt emission" 
to refer to the emission in the first peak, ~ — 13 — 1-20 s from the burst trigger (see § 12. ip . 
"flare" to refer to the peaks more than ~ 70 s after the trigger, and "afterglow" to refer to 
the smooth decay > 300 s from the trigger. 

In this paper we describe the prompt, flaring and afterglow properties of GRB 060714, 
with particular emphasis on the flares. In § [2] we discuss the observations and data analysis 
in general, while § [3] focuses on the analysis of the spectral and temporal properties of the 
flares. Finally, in § H]we show that the properties of the flares rule out an external shock (or 
long lived reverse shock) origin, and provide some (though not unequivocal) support for an 
origin common to that of the prompt emission. In §4.3 and Appendix A we exclude internal 
shocks with a small contrast in the Lorentz factor as the origin of the flares in GRB 060714, 
and point out that this mechanism faces serious problems for most Swift GRBs. 

Throughout the paper we have followed the convention F Vjt oc z/ _/3 t~ Q , where the energy 
spectral index (3 is related to the photon index T = (3 + 1. We have adopted the standard 
values of the cosmological parameters: H = 70 km s" 1 Mpc -1 , Qm = 0.27, and ft\ = 0.73. 
The phenomenology of the burst is presented in the observer frame unless otherwise stated. 
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2. Observations and Data Analysis 
2.1. Swift-BAT 

At 15:12:00 UT, 14 July 2006, the Swift Burst Al ert Telescope (BAT ) triggered and lo 



cated on-board GRB 060714 (BAT trigger=219101: iKrimm et al.ll2006al ). Unless otherwise 
specified, times t in this paper are measured from the BAT trigger time, (UT 15:12:00.3) 
hereafter designated To- The burst was detected in the part of the BAT field of view that 
was 27% coded, meaning that it was 33.6° off-axis and only 27% of the BAT detectors were 
illuminated by the source. The spacecraft began to slew to the source location at T + 34.9 s 
and was settled at the source location at T + 88.1 s. 

The BAT data for GRB 060714 between T - 240 s and T + 962 s were collected in 



event mode with 100 fis time resolution and ~6 keV energy resolution (IKrimm et al.ll2006bl ). 
The data were processed using standard Swift-BAT analysis tools and the spectra were fit 
using xspec 11.3. Each BAT event was mask-tagged using batmaskwtevt with the best 
fit source position. Mask-tagging is a technique in which each event is weighted by a factor 
representing the fractional exposure to the source through the BAT coded aperture. A weight 
of +1 corresponds to a fully open detector and a weight of -1 to a fully blocked detector. 
Flux from the background and other sources averages to zero with this method. All of the 
BAT GRB light curves shown in Figures [TJ [2J and [7] have been background subtracted by 
this method. This method is effective even when the spacecraft is moving since complete 
aspect information is available during the maneuver. 

The mask- weighting is also applied to produce weighted, background subtracted counts 
spectra using the tool batbinevt. Since the response matrix depends on the position of the 
source in the BAT field of view, separate matrices are derived for before the slew, after the 
slew and for individual segments of the light curve during the slew. 

The mask- weighted lightcurves seen in Figured] show an initial triangular-shaped (rising 
and falling power laws) peak starting at T — 15 s, peaking at T , and ending at T + 55 s. 
Given the size of the statistical error bars in the time before the slew, none of the fluctuations 
seen in Figure [TJ before the slew are statistically significant. The initial prompt emission was 
followed by two strong flares, the first starting at T + 70.4 s with a duration of ~ 10.5 s, 
and the second starting at To + 87.6 s with a duration of ~ 12 s. Finally there was a much 
weaker third flare, starting at To + 108.7 s with a duration of ~ 15.3 s. Taking into account 
the prompt emission plus the flares, we derive T 90 = 115 ± 5 s (estimated error including 
systematics). 



The time-averaged spectrum from T — 13.4 s to T + 18.0 s is best fit by a simple 
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power-law model. The power law photon index of the time-averaged spectrum is T = 
-dlogN ph /dlogE ph = 1.61 ± 0.13 (x 2 = 43.2 for 59 d.o.f.). The fiuence in the 15-150 keV 
band is (1.22 ± 0.11) x 10~ 6 erg cm -2 . The 1 s peak photon flux measured from T +75.42 
s (during the first flare) in the 15-150 keV band is 1.4 ±0.1 ph cm -2 s _1 . The prompt 
component does show spectral evolution, as evidenced by the increasing power law index: 
(T -13.4 to +2.1 s): 1.47±0.19; (T + 2.1 to +18.0 s): 1.61 + 0.17; (T + 18.0 to +70.2 s): 
2.29 ± 0.49. All the quoted errors are at the 90% confidence level. 

We attempted to fit a model consisting of a power law with an exponential cut-off to 
the prompt emission. Such a fit did not constr ain Speak, the peak of the vF{y) spectrum. 
However, it may b e possible to u se the results of IZhang et al.l (120071 ) to estimate S peak for the 



prompt emission. IZhang et al.l (120071 ) have shown that due to the relatively narrow energy 



band of BAT, it is often difficult to constrain Speak, even when Speak is within the BAT 
energy range. These authors have found that the power-law photon index T of a simple 
power law fit and S peak are well correlated with a relationship: 

logS peafc = (2.76 ± 0.07) - (3.61 ± 0.26) logT , (1) 



under the assumption that S pea k is within the BAT energy range. ISakamoto et al.l (120071 1 
have reached a similar conclusion and consistent result by combining simulations and a study 
of bursts for which Speak has been determined. We can use Equation 1 and the measured 
T = 1.61, to find that Speak = lOS.ll^s keV. This value of S pea k is consistent with the 
majority of long GRBs detected by BAT. However, another possibility is that S peak is above 
the BAT energy range (> 150 keV ), and 1.61 is i nstea d simply the low energy power law 
index a. According to the study of iKaneko et al.l (120061 ). a value of a = 1.61 is within the 



range of a values found for BATSE bursts, although only 22 of the 350 bursts in the BATSE 
sample (6.3%) have a > 1.61. Since we cannot exclude this second possibility, we will be 
conservative and quote only a lower limit, S pea k > 77.6 keV (taking the lower limit on S pea k 
as derived from Equation 1). As discussed in §EJ S peak is well constrained for the later flares. 



2.2. Swift-XRT 

The spacecraft slewed immediately to the BAT location of GRB 060714 and the Swift X- 
Ray Telescope (XRT) began observing the burst at 15:13:39 UT (T + ~ 99 s), first in Image 
mode and then Photodiode (PD) mode. The Image mode contains no spectral or timing 
information and the PD mode data is severely corrupted by the presence of two hot columns. 
Thus the fi rst us able data is in Windowed Timing (WT) mode beginning at T + 107 s (see 



Hill et al.l (120041 ) for a description of XRT readout modes). Starting from Tq+249 s all 



observations were carried out in Photon Counting (PC) mode. 
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The XRT data were processed with the xrtdas software (v. 1.7.1) developed at the 
ASI Science Data Center and included in the HEAsoft package (v. 6.0.4). Event files were 
calibrated and cleaned with standard filtering criteria with the xrtpipeline task using the 
latest calibration files available in the Swift CALDB distributed by HEASARC. Event lists 
were selected in the 0.3-10.0 keV energy band and grades 0-12 for PC mode data and grades 



0-2 for WT data were used in the analysis (see lBurrows et al.l (j2005al ) for a definition of XRT 
event grades). 

The XRT PC image of the field clearly showed a bright fading X-ray object in the field 



( IPerri et al.ll2006l ). The coordinates of the burst were determined by the XRT to be (J2000): 
RA:15 h ll m 26!5 (227?8604), Dec: -6°33'59"3 (-6?5665) with a 90% confidence error circle 
radius of 3.8 arcsec. 

Events for temporal and spectral analysis of WT mode data were selected using a 40- 
pixel wide rectangular region centered on the afterglow. Background events were extracted 
from a nearby source-free rectangular region of 40 pixel width. Data in PC mode during 
the first Swift orbit (from T + 249 s to T + 1610 s) were significantly affected by pile-up. 
By comparing the o bserved Point Spread Function (PSF) profile with the analytical model 



( IMoretti et al.ll2005l ). we removed pile-up effects by excluding events within a 2 pixel radius 
circle centered on the afterglow position and using an outer radius of 20 pixels. From the 
second orbit, the afterglow count rate was below the XRT pile- up limit and events were 
extracted using a 10-pixel radius circle. The background for PC mode was estimated from a 
nearby source-free circular region of 50-pixel radius. Source count rates for temporal analysis 
were corrected for the fraction of PSF falling outside the event extraction regions and for 
pixels partially exposed. Ancillary response files for the spectral analysis were generated 
with the xrtmkarf task applying corrections for the PSF losses and pixel exposures. The 
latest response matrices (v. 008) available in the Swift CALDB were used and source spectra 
were binned to ensure a minimum of 20 counts per bin. 



As discussed in detail in §[3], the XRT light curve (IPerri et al.ll2006l ) shown in Figures [21 
|3l and [7] displays three flares during the first orbit peaking at about T + 115, 140 and 180 
s after the BAT trigger. There is a steep decay after the end of the third XRT flare, with 
temporal power-law index a% = 2.14 ± 0.13, followed by a much shallower decay starting 
at tbrcak.i ~ 330 s, which is probably the afterglow emission. The afterglow decay from 
To + 324 s to To + 1.2 Ms can be fit with a broken power-law with an initial decay slope of 
a 2 = 0.24 ± 0.03, a break at thr^v,?. = 3.21H ks, and a post-break slope of a 3 = 1.22 ± 0.03 
(we use here the notations of Nousek et al.lboOa ). 



A power-law fit to the 0.3-10 keV spectrum from T + 107 s to T + 248 s (WT mode) 
gives a photon index of T = 2.05 ± 0.06 and a column density of (2.26 ± 0.20) x 10 21 cm -2 . 
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The Galactic hydrogen column density in the direction of the burst is 6.7 x 10 20 cm -2 . An 
extrapolation backward in time of the flat afterglow component (a 2 segment in Figure [3]) 
tells us that the contribution of the underlying afterglow during this period is negligible 
and can safely be ignored in the spectral analysis. We also fit an absorbed single power-law 
model to the XRT 0.3-10 keV spectrum in PC mode (from To + 249 s to To + 1610 ). Here 
we found a photon index of T = 2.2 ± 0.2 and a column density of (1.7 ± 0.5) x 10 21 cm~ 2 . 
At later times, from T + 6096 s to T + 45416 s, the X-ray spectrum was well described by 
a single power-law model with photon index T = 2.4 ± 0.4 and an absorbing column density 
of (1.9 ±0.8) x 10 21 cm" 2 . 



2.3. Swift-UVOT 



The Swift Ultra Violet /Optical Telescope (UVOT; iRoming et al.ll2005al ) began observing the 
field of GRB 060714 at T + 90 s in the settling mode and obtaine d the first detection in th e 
White filter (160-650 nm) in the exposure starting at T + 108 s (IBoyd fc Marshallll2006l ) . 
The position of the afterglow measured in the image from the initial exposure with the white 
filter is RA:15 h ll m 26!444 (227?8602), Dec: -e^'SS^ (-6?5662) (J2000). The uncertainty 
of this position is likely to be dominated by systematic errors, which we estimate to be ~ 0.5" 
(90% confidence radius) based o n residuals when m atching; UVOT sources in the image to 
stars in the USNO-B1.0 catalog (IMonet et al.l 120031 1 . UVOT used the standard sequence of 
exposures for observing gamma-ray bursts. The sequence cycles through all seven lenticular 
filters with increasing exposure times as the time from the trigger increases. The afterglow 
was strongly detected in the White and V filters and weakly detected in the B filter. Figure H] 
shows the UVOT detections and upper limits in the White, V and B bands, using UVOT 
C ALDB versi on 20061116, and correct ing for galactic extin ction using the extinction curve 
of |Pei I (119921 ) and the reddening from ISchlegel et al.l (119981 ) . The count rates in the White 
filter were converted to equivalent V magnitudes using the ratio of the average count rates 
in White and V seen in the multiple exposures between To+700 s. and To+1582 s. Upper 
limits (2a) in the U and Wl bands (omitted for clarity from Figure H]) are U > 18.6 (667 s 
< T < 835 s ), and Wl > 19.3 (643 s < T n < 811 s ). This lack o f detection is consistent 
with the reported burst redshift of z = 2.71 (IJakobsson et al.ll2006al See § 12.41) which shifts 
the Lyman edge to 338 nm, cutting out much of the U and all of the Wl band. The Lyman 
forest is also likely to reduce the flux in the U and B filters. 

During the time period in which the X-ray light curve is showing flares followed by a 
steep decline, the optical light curve in the V band is essentially flat and at the same level as 
the afterglow. This tells us that the flaring activity does not manifest itself in the optical and 
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suggests that the afterglow component dominates the optical light curve even at early times. 
This is consi stent with what has been seen for other b ursts with intense early flaring activity 

iGuetta et al.l 120071 ). and for the unusual "late plateau" of 
20071 ). At least during the time for which we have optical data, 



(for example iRomano et a 



GRB 070110 flTroia et al.l 



this burst seems to behave in a similar way to GRB 050401 favkoff et al.ll2005l ) in which the 
optical and X-ray emission vary independently, but the prompt optical emission is consistent 
with a backward extrapolation of the later afterglow emission. 



2.4. Other observations 



GRB 060714 was well obs erved by a large numbe r of other telescopes with detections 



(Jakobsson et al. 


2006c) 


( Jakobsson et al. 


2006d) 



2006c). A 2cr upper limit of R > 25.0 was obtained at T + 6.28 days 
2006dl ). Figure H] shows the reported optical detections compared to the 
X-ray data. During an observation beginning 12 hr after the burst, with the FOcal Reducer 
and low disp ersion Spectrograph for the Very Large Telescope of the European Southern 
Observatory, Jakobsson et al.l (12006a ) acquired a spectrum showing numerous absorption 
features which correspond to a redshift of z — 2.711 ±0.001 and a neutral hydrogen column 
density of (6.3 ± 1.5) x 10 21 cmT 2 . As seen in Figure IU the R-band optical light curve can 
be well fit to a broken power-law decay with an index a# i2 = 0.23 ±0.13 before T ± ~ 10 4 s, 
and olr$ = 1.29 ± 0.10 after. These values are very similar to the decay constants for the 
X-ray afterglow, ax,2 = 0.24 ± 0.03 and ax,3 = 1-22 ± 0.03, respectively. However the break 
in the X-ray light curve occurs at least 3 ks (corresponding to a factor of ~ 2 in time) earlier 
than the break in the R-band light curve. 



3. Spectroscopy and Light Curves of the Flares 

There were a total of five flares detected in this burst as seen in Figure [2j The first two were 
seen only in the BAT: the first was during the slew and the second during times when the 
XRT was observing in Imaging and Photodiode mode and no detailed spectral or temporal 
information is available. The third one was significant in both BAT and XRT and the fourth 
and fifth were seen as strong flares in the XRT, and weakly in the BAT. 

In the observer frame the flares lasted from ~ 84 s to ~ 210 s after the onset of burst 
emission (at T — 13.4 s). Given the burst redshift of z — 2.71 (a quite typical value for 
Swift) , in the source frame the flares extend from ~ 23 s to ~ 56 s after the start of emission. 
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Although this means that the flares are relatively early as seen in the source frame, the 
timing of flares varies widely in the cosmological frame of GRBs, and there is no indication 
that absolute timing is a distinguishing characteristic of flares. Although the timing analysis 
described in § 13.11 and discussed in §H]is done in the observer frame, the conclusions are all 
based on relative timing and do not depend on the cosmological frame chosen. 



3.1. Spectroscopy 

All five flares were fit with a spectral model of a power law with exponential cut-off: F(E) = 
A(E/100 keV)~ a exp[— E(2 — a)/E peak ], where E is the photon energy, E peak is the peak 
energy of the vF{y) spectrum, a is the photon index, and A is a normalization factor. The 
first two were fit to BAT alone, the last three to BAT jointly with XRT. For the fourth and 
fifth flares inclusion of the BAT data did not significantly affect the fit. The bars at the 
bottom of Figure [2] indicate the time segments used to fit each of the flares. The results of 
the spectral fits are shown in Table [T] and in the top two panels of Figure Although there 
is no evidence for a smooth power law decay underlying the flares, it is possible that they 
overlap each other temporally. We have increased the error bars on the flux and Ei so values 
in Table [T]to account for this overlap by extrapolating the power-law fits (§ 13. 2p to each flare 
down to zero and estimating the fraction of the flux falling outside of the nominal start and 
stop times of the flare (for the upper limit), and the flux possibly due to neighboring flares 
(for the lower limit). 

In the fits to the last three flares an absorption component was included in the model. 
The column densities were found to be (units 10 21 cm -2 ), respectively for flares 3-5, 1.91 ± 
0.43, 1.86 ± 0.36, 1.67 ± 0.22, consistent within errors to a constant value. 

One can see in Figure [5] that the peak energies of the flares decrease with time and 
the power law indices show a general softening of the spectra. Furthermore, the apparent 
linearity of the plots indicates a connection (or at least a clear trend) between the five flaring 
events. The time dependence of E pea _ k is well fit by a power law, with index —5.81 ± 0.68. 
Similarly, the time dependence of the power law index can be fit to a power law with index 
0.67 ±0.15. 

Using the redshift z = 2.71 for this burst, we extrapolate the total isotrop ic equivalent 
radiat ed energy, i? 7) j SO (in ergs) in the range of 1-10 4 keV, using the definition of lAmati et al. 



( 120021 ). For the flares, the extrapolation fixes the -E pea k and a values derived from the cut-off 
power law fits, and uses a fixed high energy index (3 = —10.0. For the prompt emission we 
derive a lower limit to -E 7i i S o from the lower limit to -E pea k- The third panel of Figure [5] shows 
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-E 7i i SO as a function of time. There is a general trend towards lower total energy output of 
successive flares. Here a fit to a temporal decay power law gives an index of —1.72 ± 0.46. 
In Figure [6] we plot E peak against -E 7i i SO for the flares and the prompt emission. This enables 
us to compa r e the episodes of GRB 0607 14 to the E peflk — E^^ n relationships found by 



Amati et al.l (120021 ) and parameterized by iGhirlanda et al.l (12004 ) The prompt emission 



and the first two flares fall on or very close to the relationships, while the last three flares fall 
well below them. This shows that -Epeak is falling with time more rapidly than the square of 
the total isotropic equivalent energy emitted in the flares, (-E^iso) 2 - 



3.2. Temporal Analysis 

In addition to the spectral properties of the flares, their temporal properties can also 
help to tell us whether or not these events arise in the external shock (i.e. are afterglow 
emission). 

In order to study the fine time structure of sub-peaks within the flares, we have derived 
a robust method to distinguish between a significant change in slope of the light curve 
(indicating the start or end of a new subpeak) and a statistical fluctuation. Using one- 
second binned light curves for both the BAT and XRT data, we applied an iterative method 
to find each significant episode (subpeak) during the times of the flares . In the first iteration, 
peaks were defined as local maxima (points higher than each of their nearest neighbor points) 
and valleys as local minima. In successive iterations we culled the peaks by requiring that 
a significant peak be at least three standard deviations above the valleys on either side. By 
this method, each point in the light curve was assigned to a specific interval: either the rise 
or fall of a peak or to the periods of slow rise (T < —13.4 s) or slow decay (18.0 s < T < 70.2 
s or T > 105.4 s for BAT; and T > 195.6 s for XRT). We then fit each interval to a power 
law: log(R) vs. log(t), where R is the count rate for either BAT or XRT (depending on 
which flare is being analyzed) and t is the time since the burst trigger. Finally we use these 
power law fits to define the actual start, apex and end of each subpeak (not restricted to light 
curve bin edges). The start of each episode (ti) is defined as the time at which the rising 
power law segment of a peak crosses the falling power law segment of either the preceding 
smooth decay or the previous peak. Similarly the peak time of each interval (£2) is the time 
at which the rising power law segment of a peak meets the falling power law of the same 
peak. The rise time (At) is thus defined as At = t% — t\, while the time associated with each 



1 Since we do no t see a jet break in the light curve, we are unable to constrain the -E pea k — E 1 (9) relation 
of IGhirlanda et alj (|20o4) . 
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rise episode is t = (t 2 + ti)/2. Thus 

A*_2(t 2 -*i) At At 



(2) 



t (t 2 + h) (h + At/2) (t 2 -At/2)' 

The conclusions drawn in § H] do not depend critically on the exact definitions of t or At. 

The episodes so defined are shown graphically in Figure [7] and their main temporal 
properties in Table [2j The temporal decay indices are derived in two ways, detailed in the 
table caption. The decay index is first calculated using the burst trigger time To as the 
reference time (a a), which is the standard way that GRB decay indexes are calculated, 
and would be appropriate for the flares if they were afterglow emission. However, since 
the resulting values of the decay index a a are very high, cor responding to a very steep 



decay which is very hard to produce by the afterglow emission (IKumar fc Panaitescul 12000 



Nakar fc Piranll2003l ). we also want to explore the possibility that the flares arise from late 
time sporadic activity of the central source. In this case the appropriate reference time, t , 
would roughly correspond to the onset of the individual flare or sub-flare whose decay rate 
we wish to quantify (cub)- Both decay indices are shown in Table (2], although values quoted 
in Figure [3] are taken as the decay index q b . Even taking the second definition, q b , the decay 
indexes of the flares are all very steep, and except for flare 4, much steeper than either the 
decay of the emission just before the first flare (a = 1.16 ± 0.47) or the decay immediately 
after the fifth flare (a = 2.14 ± 0.13). The decay slope immediately after the last flare 
is still much steeper than the afterglow beginning at T + ~ 320 s, so it is likely part of 
the prompt emission as well. We note, however, that «b does not exceed 2 + (3 (within the 
statistical uncertainty, where F„ oc v~^{t — tn)~ a ) : whic h is the steepest decay allowed by 



the 'high latitude' emission (IKumar fc Panaitescul l2000l ). and thus the decay of the flares 



and sub-flares is consistent with the expectations for late time intermittent activity of the 
central source. 



4. Discussion 



4.1. The External Shock 



The temporal properties of the flares provide strong evidence against an external shock 
(i.e. afterglow) origin for them. Figure [H] shows the fractional increase in flux, AF/F, versus 
the ratio of the rise time (At) to the peak of each flare or sub-flare and the time (t) from 
the GRB trigger, At/t (see Eq. [2]). It can clearly be seen that large amplitude variations in 
the flux, AF/F > 1, occur on very short tim e scales, At/t <C 1. This basically rules out an 
external shock origin for the flares (see, e.g., Iloka et al.ll2005l ; lNakarll2006l ; lLazzati &: Perna 
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20071 ). We also note that the flares are very different from the smooth, late time tail emission 
in the 20-100 keV band which was observed to foll ow many of the bur sts detected by the 
Burst and Transient Source Experiment (BATSE) (jConnaughtonl 120021 ) . The tail emission 
is temporally much smoother than the flares of GR B 060714, an d at a much lower level (as 
a ratio to the peak of the prompt emission). Also iGiblin et al.l (120021 ) identify a subset of 
the BATSE bursts that have high-energy decay emission consistent with forward external 
shocks. Although it is below detectability for BAT, it is possible that the tail emission is seen 
after the last flare in the ot\ = 2.14 decay segment (Figure [3]) • This decay index is consistent 
with the average t emporal decay index reported for BATSE tail emission, —2.03 ± 0.51 
jGiblin et alJl2002h . 



The major possible sources of variability in the afterglow light cu rves (i.e. in the emis- 



sion from the external shock) are: (i) a variab le external density (jWang &: Loebl 12000c 
Lazzati et al.l 120021 ; iNakar. Piran fc Granotl 120031). fii) a "patchy shell" i.e. angular inho- 
mogeneity within the outflow (IKumar fc Piranl l2000bl ; INakar. Piran fc Granotl 120031 ) , and 
(hi) "refreshed shocks" i.e. relatively slow shells that were ejected from the central source 
toward the end of the prompt emission and catch up with the after glow shock as the latter 



decelerates to a Lorent z factor slightly lower than that of the shells (IRees fc Meszaroslll998 
Kumar fc Piranl l2000a ). Such a sharp (At/t <C 1) large amplitude (AF/F > 1) rise in 



the observed flux, as we find for GR B 060714 (see Figure El) , cannot be caused by a sud- 
den increase in the external density (|Nakar &: Granotl 120061 ). In addition, a "patchy shell" 
produces At/t ~ 1 (INakar &: Oreni 120041 ) . and cannot account for the observed At/t <C I, 
since new 'bright spots' in the outflow become visible (i.e. enter the observed region of 
angle ~ I/7 around the line of sight) gradually, on the dynamical time (At ~ t). Finally, 
"refreshed shocks" produce At/t ~ 1 before the jet break time, when the rise time At is 
dominated by the angular time tg ~ R/2cy 2 . After the jet break time, the angular time is 
tg ~ R9 2 /2 (where 9j is the half-opening angle of the jet), and it soon becomes smaller than 
the radial time t r ~ R/lOcy 2 (since 9j remains close to its initial value, 9 , as long as the jet 



i s relativistic; lGranotll2006f). s o that the rise time At becomes dominated by the radial time 
(IGranot. Nakar &: Piranll2003l ). However, even in this case At/t > 0.2, which cannot explain 
the values of At/t < 0.1 and in some cases even as low as At/t < 10 -2 , that we obtain in 
our analysis. 



4.2. Other Possible Causes for the Flares 



Now that we have effectively excluded an external shock origin for the flares, we explore 
other possible explanations. These involve sporadic late time dissipation events which are 
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a result of either (i) early time activity of the central source on a time scale comparable 
to the duration of the prompt GRB emission, or (ii) late time intermittent activity of the 
central source, which is rather directly reflected in the observed times of the flares. Both 
types of models may in principal be applicable both for the prompt GRB emission and for 
the flares, due to their roughly similar observed properties. Thus we also wish to address 
the similarities and differences in the observed temporal and spectral properties of the flares 
and the prompt emission. 

Models in the literature that can be included in the first of the two above classes include 
both emission powered by sporadic magnetic dissipati on events within the outflow, possibly 
induced by its interact ion with the external medium (ILyutikov fc Blandfordl 120021 ; iGiannios 
20061 ; iThompsonl 120061). and late time internal shocks between shells with a small relative 
velocity (e.g., iBarraud et al.l 12005c iBurrows et al.l 120071 ). In the former models the temporal 
and spectral properties of the emission have not yet been worked out in detail, so direct 
comparison to observations is not possible at this stage. We now briefly address the latter 
model. 



4.3. Internal Shocks with Small Contrast in the Lorentz Factor 

In this scenario the difference in the Lorentz factors of the colliding shells, A7, is much 
smaller than their typical Lorentz factor 7: A7 <C 7. For convenience, we provide most of 
the relevant results here, while a detailed derivation of these results as well as some relevant 
(but somewhat more technical) discussion is provided in Appendix [A] In this picture, later 
collisions correspond to a smaller A7/7, occuring at an observed time tfl are and with a 
duration At flare which satisfy 

1 1 

^flare ~ £ e j + A£ e j , At flare ~ ^— At e j ~ t flare — t e j , (3) 

where the first and second shells are ejected at times t e j and t e j + At e j with Lorentz factors 7 
and 7 + A7, respectively. This makes it hard to account for the short time scale variability 
At/t <C 1, unless the colliding shells were ejected from the source at a time t e j after the GRB 
trigger that is much closer to the observed time of the flare tfi are than to the GRB trigger 
t = (since At/t = Atfl are /tfl arc ~ 1 — t e j/tfl ar e)- However, this corresponds to models of class 
(ii) above (where the central source is active a t late times, close to the time of the observed 



flares), rather than class (i). lLazzati fc Pernal (120071 ) come to a similar conclusion. 



In the latter case, which corresponds to class (ii), one can in principal account for the 
temporal properties of the flares. We also wish to examine whether their spectral properties 
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and energetics can naturally be reproduced. In this picture the shells collide at a radius 

#is^^7 2 cAt ej , (4) 

(the subscript 'IS' is for internal shocks) which is larger by a factor of ~ 7/A7 ^> 1 (for the 
same average 7) compared to internal shocks with a reasonably large contrast in the Lorentz 
factor, A7 > 7. For a reasonably large contrast in the Lorentz factor (A7/7 > 1) the peak 
of the vF{y) spectrum, E peak , typically corresponds to hv m where v m is the synchrotron 
frequency of the relativistic electrons with the minimal random Lorentz factor in the power 
law distribution of energies. However, for low contrast internal shocks (A7/7 <C 1), 

u m oc (^) " iglg oc (^-) 6 L£ 7-^.,)-' , (5) 

where L iso is the isotropic equivalent kinetic luminosity of the outflow, under the standard 
assumptions that the fractions of the internal energy behind the shock in the relativistic 
electrons (e e ) and in the magnetic field (e#) are constant. That is, v m decreases very rapidly 
for a small contrast A7/7. On the other hand, the cooling break frequency v c increases as 
A7/7 decreases, 

,„ oc (£2) " L^R m y° oc (^) L^A U . (6) 

Therefore, for A7/7 <C 1, v c > v m an d there is slow cooling, i.e. -E pea k = hv c rather than 
hu m . 

For a reasonable value of the magnetic field that is advected with the outflow from the 
central source, such a field would dominate over a sub-equipartition shock-generated field in 
the shocked regions of the colliding shells (see Appendix [A] for details). In that case u m /u c 
scales "only" as (A7/7) 6 , instead of (A7/7) 10 for a field that is a constant fraction of the 
equipartition value, that was assumed above. Since this is still a very high power of A7/7, 
this does not change the main conclusions. 

The very strong dependence of v m jv c on A7/7 further decreases the radiative efficiency, 
which is already low for A7/7 <^ 1, since the fraction of the total energy of the colliding shells 
that is converted into internal energy (out of which the fraction that goes into relativistic 
electrons, e e , may be mostly radiated away for v c < u m but not for v c > u m ) is given by 

where x is the ratio of the rest masses of the two colliding shells (for a fixed contrast of 
the Lorentz factor, A7/7, the efficiency e is maximal for equal mass shells, x — 1). For 
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GRB 060714 the combined -E 7i i S o of the five flares is comparable to that of the prompt emis- 
sion, and that of individual flares is at most one order of magnitude smaller. This suggests 
that the kinetic energy that remains in the colliding shells that produce the flares in this 
scenario is much larger than that of the original shells that produced the prompt emission. 
This would result in very significant episodic energy injection into the afterglow shock, i.e. 
"refreshed shocks" which are inconsistent with the temporal properties of GRB 060714 (and 
may prove problematic for most Swift GRBs with X-ray flares). Furthermore, since the total 
-Ey.iso of the flares is ~ 10 52 5 erg, such an extremely inefficient emission would imply a huge 
remaining isotropic equivalent kinetic energy, in excess of 10 55 ' 5 erg for A7/7 < 0.1, which 
would imply a very large total kinetic energy (> 10 53 erg), even when corrected for the 
fractional solid angle occupied by the jet (/& > 10 -2 5 ), given that there is no sign of a jet 
break at least until ~ 10 6 s. 



4.4. Internal Shocks with Reasonably Large Contrast in the Lorentz Factor 



We now examine the more popular version of the internal shocks model, where the 
difference in the Lorentz factor of the colliding shells, A7, is of the order of their average 
Lorentz factor 7, i.e. A7 > 7). In this case the observed time of the flares reflects the time in 
which the colliding relativis tic shells were ejected from the source (jKobayashi. Piran fc Sari 



1997 



N 



akar &: Piran 



source (IBurrows et al 



200~2L thus requiring intermittent late time activity of the central 



2005bJ). Such a sporadic late time activity of the central source may 



arise (e.g., iPerna et al.ll2006l ) by infall of material into the central black hole from instabil- 
ities in an accretion disk, or by late t ime fall-back of m aterial in the collapsar model (e.g., 
MacFadyen. Woosley fc Hegerll200ll ). iFan fc Weil (120051 ) show that this model is consistent 
with X-ray flares in a number of different bursts. 



In this scenario the radiative efficiency may be reasonably high (but still typically 



< 



10%), and -Epeak is usually given by hv m . All of the equations of § 14.31 still remain valid, 
up to a factor of order unity, under the substitution A7/7 = 1. Thus Ai flarc ~ At e j and 
E peak oc L 1 1 / 2 7" 2 (A)f: e j)" 1 cx L 1 1 / 2 7" 2 (At flare )" 1 . Figure [9] shows that the rise time At ~ At flarc 
increases with time (by about two orders of magnitude), while At/t increases by a somewhat 
smaller factor (but still a factor of ~ 30 over less than a factor of 3 in time). This alone 
can roughly account for the decrease in -E^k (see upper panel of Figure E]), which decreases 
by about two orders of magnitude, similar to the increase in At flare . Since for each flare 
L iso ~ i? 7] i SO /Atflare, and -E 7i i SO of the individual flares decreased by about one order of 
magnitude, we can infer that L; so decreased by about three orders of magnitude, which 
requires 7 to decrease by a factor of ~ 10 3//4 ~ 5 — 6. Since Ris ~ 7 2 cAt e j ~ 7 2 cAtfl are , this 
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implies a modest increase in R\s by a factor of 
while Atfl are increased by a factor of ~ 100). 



3 (since 7 2 decreased by a factor of ~ 30 



4.5. Possible Relation to Models for the Late Time Central Source Activity 

The sharpness of the peak features can be seen in Figure in which the rise time At and 
the ratio At/t are plotted against the time since the trigger. There is a tendency for peaks 
to become less sharp (increasing At) with time. Th is trend is consistent with the viscous 
disk evolution model discussed by lPerna et al.l (120061 ). in which both the accretion time scale 
At and the arri yal time t depend on t he radial location of the accreting material within the 
accretion disk. iKocevski et al.l (120071 ) show that the correlation of broader pulse durations 
with time holds in general even among flares from different GRBs, supporting the idea that 
later flares are coming from larger radius shells. The lower panel in Figure [9] shows that even 
for the late flares with longer rise times, that At/t remains very small, < 0.1 for all peaks, 
although At/t tends to increase with time (it reaches values < 10~ 2 at early times). 

The distribution of At/t yalues for GRB 60714 can be compared t o the statistical 



samples of flares presented by IChincarini at al.l (120071 ) and iBurrows et al 



(2007). We see 



that m ost of the points in our At/t distribution are consistent with the IChincarini at al. 
(120071 ) sample, while the three points with At/t < 0.01 lie outside the distribution. This 
is explained by t he dif ference between our definition of At (rise time) and that used by 
Chincarini at al.1 (120071 ) (a of a Gaussian fit to the peak). Although these measures are of 
the same order of magnitude, our definition leads to very small values of At/t for those 
peaks with fast r i ses an d slower decays, and in general a broader distribution than seen by 
Chincarini at al.l (120071 ) . Given the asymmetric shape of many of the peaks, the rise time 
may be a true r meas ure of the variability, since the width can be skewed by a slow decay. 
Burrows et al.l (120071 ) show At rise /t, where At rise is defined as the time between when the 
power law crosses the underlying afterglow and the peak time, a definition which would 
lead to somewhat larger values of At r i se /t than ours. Their distribution, however, has a 
median value ~ 10 times larger than ours, indicating that the time st ructure in GRB 60714 
is unusually sharp compared to most other flares. The analysis of IKocevski et al.1 (120071 ) 
also shows that GRB 060714 has shorter flare durations and rise times than most bursts. 
Alternatively, it could be that we find significant structure on smaller time scales than in 
other works because we are analyzing brighter flares with better photon statistics. We are 
also looking at sub-structure in the flares, as is often done with the prompt emission, in 
which significant temporal structure is usually found down to the smallest time scale that 
can be measured (which is limited either by photon statistics or the temporal resolution of 
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the instrument) 



The successive flares in GRB 060714 are each spectrally softer than the preceding flare, 
as seen in the top panel of Figure The spectral softening of the flares is supportive of 
models in which m ultiple masses accrete onto a central black hole at successively later times. 
King et all (120051 ) discuss a model where the collapse of a rapidly rotating stellar core leads to 
fragmentation and the formation of multiple compact objects which accrete onto the central 



black hole on time scales which dep end on the orbital radii and fragment masses. iKing et al. 



( 120051 ) and iBurrows et all (j2005bl ) suggest that successively later accretion events occur in 
cleaner, lower density environments, since earlier jets would have excavated channels through 
the progenitor s tar. This suggests that later outflows would have lower baryon loading and 
hence higher 7. King et al.l (120051 ) also suggest that tidal effects may smooth later accretion 
events, lengthening At e j. Both increasing 7 and increasing At e j would lead to a reduction in 
-Speak as a function of time (see § 14. 4p . Note that -E pea k for the prompt emission is slightly 
higher than -E pea k for the first flare, which would be expected since the initial collapse should 
have the lowest bulk Lorentz factor. However, as demonstrated in § |4.4[ while Figure [9] 
indeed shows that At ~ Atfl are increases with time, we find that 7 should have decreased 
in time (by a factor of ~ 5 — 6) rather than increased with time (for the standard internal 
shocks model) as expected in the above model. 

As can be seen in the third panel of Figure -E 7i j SO decreases with time roughly 
as ~ t -1 ' 7 , and L iso decreases roughly as t~ 3 . This is a steeper decline than the ex- 
pected late-time accretion rate due to fallback in the collapsar model, M acc oc t -5 / 3 (e.g., 



MacFadyen. Woosley fc Hegerll200ll ). and in fact drops roughly as Mf cc , suggestive of a rel- 



ativists outflow powered by neutrino-anti neutrino annihilation, where L iso oc L v oc M^ cc 
(this is, of course, highly speculative at this stage). Furthermore, such a steep decay of L iso 
is muc h steeper than oc t 9 with q > — 1 that is required in order to account for the flat decay 
phase (INousek et al.l 120061 ) . Although the latter requires a temporally smooth outflow, and 
the flares require intermittent source activity, it is hard to see how the overall temporal decay 
of such a smooth component and a variable component would be that different. 



4.6. Similarities and Differences between the Flares and the Prompt Emission 

A striking feature of the spectroscopic results is that not only is E peak declining with 
each flare, but the low energy spectral index a is increasing (second panel of Figure [5]), 
meaning that the flares are becoming successively softer with time. This is consistent with 
the third panel of Figure [5] which shows the the total isotropic energy is decreasing with 
time. Hard-to-soft spectral evolution is commonly seen in the prompt emission of GRBs 
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( IBand &: Ford! Il998l ; iNemiroff et all Il994l ). and the bottom panel of Figure [5] shows this 
trend very clearly for the period before the start of the first flare (first three points). For 
the flares, not only is there a softening from flare to flare, but within each of the first four 
flares there is a tendency for the spectrum to soften as well. That this well-established 
feature of the prompt emission of GRBs is see n within and across the fl ares favors similar 
origins for the flares and the prompt emission. iButler fc Kocevskil (120071 ) have shown for a 
sample of 27 SW/%-detected GRB flares, that flares exhibit significant hardness-intensity and 
hardness-fluence correlations which match closely the correlations observed for GRBs. They 
attribute this hardness evolution to an evolving E pea k- The hardness ratio H = S (50-100 
keV)/S(25-50 keV) is close to unity for the prompt emission and the first two flares, so all 
three of these episode s fall well within th e distribution of long bursts in a GRB hardness- 
duration diagram (see lRoming et all 120061 for a recent BAT hardness- duration distribution 
plot). If the overall spectral shape of the flares remains roughly the same, we would expect 
the hardness ratio H to fall dramatically as E pea k shifts downward through the BAT energy 
range (top panel, Figure [5]), since the BAT spectral index would shift from the super- E pea k 
value of ~ 1 to the sub-E peafc value of ~ 2.5. This is indeed what we see in the bottom panel 
of Figure [5J supporting a common origin for the flares. 

Another interesting effect is seen in Figure [6j Here we see that the prompt emission 
and the first two flares sh ow a relation s hip b etween .Epeak and E 7iiso quite consistent with 
the relationships found by lAmati et al.l (120021 ) .This supports a similar mechanism for these 
flares and for the prompt emission, since it suggests that these flares are spectrally and 
energetically like weak, low -E pca k bursts. The last three flares fall well below the relations 
since for the succession of flares Epeak is falling more quickly than E lj - mo , with these flares 
remaining relatively energetic (E lj j so > 10 51 erg, even as -E pea k decreases to < 1 keV). This 
argues against a similar origin for these two flares and the prompt emission (especially since 
almost all known outliers to the Amati relation are on the other side of it, i.e. with a higher 
-Epeak for their E 7iso than expected by this relation, while here the situation is the opposite), 
although there are still possible ways around this. Alternatively, one could certainly argue 
that such a wide dispersion within a single burst over a wide range of energies should call 
for a re-evaluation of the Amati relationship. Although the E„ P ^ — E 7 iso relations have 



been shown to hold for many X-ray flashes (e.g.. ISakamoto et al.ll2006l ). it is much less clear 



whether the relations should hold for late-time flares or individual peaks within bursts. 

The time structure of the flares is quite similar to that of the prompt emission. All of 
the flares show structure on time scales of ~ 10 _1 — 10 s within an overall envelope of 
duration ~ 10 — 50 s. The prompt emission falls within a ~ 30 s envelope, and although 
the sub-peaks in the prompt emission are not formally significant, there is a suggestion of 
structure on shorter time scales, especially in the 50-100 keV band (middle panel of Figured]). 
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Certainly short time-scale variations in the prompt emission are a common feature of GRBs. 



5. Conclusions 

The temporal properties of the series of five flares in the long gamma-ray burst GRB 060714 
provide strong evidence against an external shock (or a long lived reverse shock) origin for 
these X-ray flares, and are consistent with sporadic late time activity of the central source. 
The strongest argument against an external shock origin of the flares is that large amplitude 
variations in flux (AF > F) occur on very short time scales (At/t <C 1, where initially 
At/t < 10~ 2 , while for all flares or sub-flares At/t < 0.1; see Figs. [S]and[HD- As discussed 
in §4.1[ this cannot be reasonably accounted for by emission from the forward shock, and 
similar considerations apply also for a long lived reverse shock. 

In the context of the internal shocks model, we show that the temporal properties of 
the flares (At/t <C 1) exclude scenarios where the colliding shells are ejected from the source 
well before the observed times of the flares (see § 14. 3j) . and require instead intermittent late 
time activity of the central source. Even in the latter picture, an internal shocks model with 
small contrast in the Lorentz factor is inconsistent with the observations (see § 14. 3j) . for the 
following reason. The large -E 7i i SO of the five flares, together with the small efficiency of such a 
model in converting the bulk kinetic energy into the observed radiation, requires that a large 
amount of kinetic energy remain in the shells, a condition that would unavoidably produce 
very prominent "refreshed shocks." Such episodic energy injections into the afterglow shock 
are, however, inconsistent with the temporal properties of the smooth X-ray (and optical) 
light curve that follows the flares. 

An internal shocks model involving sporadic late-time activity of the central source and a 
reasonably large contrast in the Lorentz factors of the colliding shells (§ I4.4[) is more consistent 
with the data. This type of model can avoid prominent "refreshed shocks," in addition to 
being able to account for both the observed temporal properties of the flares and the decrease 
with time of their -E pea k and E 7j i B0 . Although the data from GRB 060714 are not sufficient for 
us to conclusively discriminate between particular models of late-time intermittent activity 
of the central source ( §4.51) . they can start testing and perhaps eve n eventually discri minate 



between the different models. For example, the viscous disk model (IPerna et al.ll2006l ) can in 
principal explain the trend of peaks be coming less sharp with increasing time . The simplest 



expectations of fragmentation models (IKing et al.ll2005l ; iBurrows et al.ll2005bl ). however, are 



axp 

m 



inconsistent with the decrease in the Lorentz factor 7 that follows from the time decay of 
i?pe a k and -E 7i i SO in the internal shocks model. Moreover, the decline in E^^n is steeper than 



what would be expected from fallback in the collapsar model (IMacFadyen. Woosley Heger 
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200ll ). but suggestive of a relativistic outflow powered by neutrino-anti neutrino annihilation. 



Some evidence in favor of a common origin for the flares and the prompt emission is 
provided by the spectral and temporal similarities of these components. As in the prompt 
emission, there is a clear spectral softening as the flares progress, as well as within most 
of the flares. Furthermore, both the prompt emission and the flares show time structure 
on the scale of ~ 1 s within an overall envelope of duration ~ 10 s. On the other hand, 
an arguments agai nst a common origin is the inconsistency of the final three flares with the 
Amati et al.l (120021 ) relation, although it is not well-established that this relation should hold 
for late-time flares. 

After the fifth and final flare, GRB 060714 showed typical afterglow time profiles con- 
sisting of three separate power law segments (Figure 13} • This part of the light curve is a 



nice example of the "canonical" afterglo w light curve found by Swift (INousek et al.l 12006 



Panaitescu et al.ll2006l ; IZhang et al.ll2006l ). This burst was well observed at late times in the 
optical (Figure H]) and the optical light curve shows both similarities and differences with 
respect to the X-ray light curve. The first optical observations showed emission at a typical 
value for an optical afterglow in the Swift era (V = 18.6 at ~ 200 s), and the optical magni- 
tude remained roughly constant out to at least To + 7700 s, with no indication of an optical 
counterpart to any of the flares. Similar to the X-ray light curve, the R-band light curve 
also had a break to a steeper temporal decay, although the optical break occurred about a 
factor of ~ 2 in time later than the X-ray break. Both before and after the break, the X-ray 
and R-band power law temporal decay indices were quite similar. This suggests a similar 
origin for the late-time X-ray and optical emission — most likely afterglow emission from 
the forward shock, probably with some contribution from a long lived reverse shock before 
and around the time of the break, which might potentially account for the difference in the 
breaks times between the optical and the X rays. 

The subject of X-ray flares in GRBs is an important one and a number of recent papers 
on Swift bursts have discussed flares, leading to a growing consen sus that flares reflec t late 
time activity of the central s ource. Some import ant examples are lFalcone et al.l (12006) (the 
gia nt flare of GRB 050 502B) jGuetta et all (bo07h (GRB 050713A). and lBurrows et alifhoOTt ) 
and lChincarini at al.l (120071 ) on the analysis of a statistical sample of Swift flares. This paper 
adds important new observational results and conclusions to the existing body of literature 
on flares. GRB 060714 contains an unusually large number of well monitored flares in a 
single event, and our interpretation of the timing and energetics of these flares not only 
provides strong evidence against an external shock or a long lived reverse shock model, but 
allows us to favor the model of late-time activity of the central source and a large contrast in 
the Lorentz factor of the colliding shells over an internal shocks model with a small contrast 



in the Lorentz factors and refreshed shocks. This event provides the best evidence to date of 
a continuous and gradual transition in the spectral and temporal properties of the prompt 
emission and flares, which is very suggestive of a common origin. 
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A. Derivation of the Equations in § 14.31 

Here we provide a derivation of the equations for internal shocks with a small contrast 
in the Lorentz factor, that appear in § 14.31 Consider two shells, where the (back edge of the) 
first shell (subscript '1') is ejected from the source at (lab-frame) time t\ = £ e j with Lorentz 
factor 71 = 7 and rest mass mj, while the (front edge of the) second shell (subscript '2') 
is ejected at t 2 = t e - s + At e j with a Lorentz factor 72 = 7 + A7 and rest mass m 2 . We are 
interested in the limit of low contrast in the Lorentz factors of the two shells, A7/7 <C 1, 
and provide below approximate expressions which are valid in that limit. The two shells 
collide at a lab frame time tis and a radius Ris which satisfy Ris = -Ri(tis) — -^(^is) where 
Ri(t) = Pic(t — t e j) and i?2(£) = @2c(t — t e j — At e j) so that t\s = £ e j + AjAt B j/ (P2 ~ Pi) an d 

P2-P1 1 - (71/72) 2 A7 

where the "=" is valid for 71, 72 3> 1, while the requires also A7/7 <C 1. 

The observed time corresponding to the onset of the resulting spike in the light curve 
(or "flare") is the arrival time of a photon emitted along the line of sight at Ris and t\s, 
relative to a photon emitted at R = and t — 0, which is given by 

W.onset - tlS " — " + ^ _ A " ^ + - _ (7i/72)2 ~ t eJ + ^ -y ■ {A2) 

The exact observed time at which the flare peaks, tfl are , would be somewhat later, typically 
by about the shell shock crossing time. The latter depends on the details of the shell (width, 
Lorentz factor and mass density distribution within the shell, etc.) but is generally expected 
to be of the order of t n are,onset — £ e j • Furthermore, the typical width of the spike in the light 
curve, Atfl are , is of the order of the angular time (which is typically also of the order of the 
radial time or shock crossing time), and since Ai?ig ~ Ris we have 



Ris 7 

Atfl are ~ ~ ~/\^^ C i ~ _ ^ ~ ^ (^flare,onset ~ *cj) • (A3) 



In order for A7 to be meaningful, we assume that it is larger than the spread in the Lorentz 
factor within each shell, so that the shells do not spread significantly before Ris and their 
width (in the lab frame) is of order cAt e j. 

The values of the spectral break frequencies depend on the physical conditions within 
the shocked shells. The isotropic equivalent kinetic luminosity of the outflow is given by 
L iso « 47ri? 2 7 2 p'c 3 while the total luminosity is (1 + a)L iso where a = (B') 2 /Airp'c 2 is the 
ratio of electromagnetic to kinetic energies of the outflow, and primed quantities are measured 
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in the comoving frame. The relative velocity between the two shells is 



Al= A^L^2|^4«^2« 1 . (A4) 
P2P1 72 + 7i 7 



Thus, the shocks going into the two colliding shells are Newtonian for A7/7 <C 1, and the 
compression ratio is p' ps / p' « 4, where the subscript "ps" is for post-shock. For equal density 
shells, the relative velocity of the upstream and downstream fluids is /3 u d ~ P21/2 ~ A7/27, 
and the internal energy per unit rest-energy is El e' ps / 'p' ps c 2 ~ (3 2 d /2 m (A7 /^y) 2 /8. Thus, the 
minimal random Lorentz factor of the power-law distribution of relativistic electrons scales 
as 

7 m oc e e — oc e e I 1 , (A5) 

Pps V 7 / 

where e e is the fraction of the post-shock internal energy that goes into such a population of 
relativistic electrons. The internal energy in the shocked regions scales as 

~ s {-) » c « {-) 7w « S 7(A(-f ■ (A6) 

where we evaluate the values of the relevant quantities at R w i?is, using Eq. IAU 

If the magnetic field holds a constant fraction, e#, of the internal energy behind the 
shock, as is often assumed for a shock-generated magnetic field, then such an equipartition 
field □ would scale as 

« ( £Be y 1/2 oc 4 /2 (-^) 2 LVV 3 (A iei) - . (A7) 

However, some magnetic field is expected to be advected with the outflow from the central 
source. At large distances form the source such a field is expected to be primarily in the 
tangential direction (normal to the radial direction) so that it would be amplified in the 
shock by the compression ratio, -B^dvps ~ ^-^adv where (-B^dv) 2 = 47rp'c 2 cr. Thus, for equal 
mass shells the ratio of the magnetic energy density associated with this field to the internal 
energy density in the shocked region is (-B^ dv ps ) 2 /87rep S ~ 64a(A7/7)~ 2 , so that the shock 



2 In this case, if the mass of the shells is also the same, the two shocks finish crossing the two shells 
together, and this is also the fraction e of the total energy that is converted into internal energy. As shown 
below, for a fixed Lorentz factor contrast A7/7, e is maximal for equal mass shells. 

3 We use the term "equipartition field" for simplicity, even though strictly speaking it is holds a constant 

1 /2 

fraction (e B , generally smaller than unity) of the equipartition value. 
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compressed advected field would typically exceed an equipartition shock generated field for 
reasonable values of a, i.e. a > 10 _3 (A7/7) 2 (e B /0.1). In this case B' ps « -B^ dvps , where 

^dv, ps « Scv^V oc L^V^Atej)- 1 • (A8) 

The synchrotron frequency of the electrons with the minimal random Lorentz factor, 7 m , 
scales as 

6 



eB> ^ [ (¥) ^V^)" 1 (equipartition) 
Vm « 7 pL>™ K J V ^ (A9) 

" ' '"' ' ! rV2 (^r) e ^i/ 2 7 ^(At ej )^ (advected field) . 



The random Lorentz factor of electrons that cool on the dynamical time (shell crossing time), 
t' ~ -Ris/7, scales as 



6irm e c j B V 7 



^T 1 ^ 1 (^) " ^o7 5 At ej (equipartition) 
1 + y)- 1 ^ 1 L^V A i e j (advected field) 



(A10) 



where Y is the Compton y-parameter. 



The factor of (1 + Y) can be safely dropped from the expression for the advected field, 
since in this case Y <C 1 (it is included her e for completeness ) . This can be seen as follows. In 
general Y(l+Y) ~ /3 2)S he r ade e / £b (see, e.g. lSari fc Esirul200ll ) where j3 2tS h is the velocity of the 
downstream medium relative to the shock front, which in our case is (3 2 ,sh ~ P21 ~ A7/7 <C 1, 
and e ra( j ~ min[l, (v m j V c ) ^~ 2 )' 2 ] is the fraction of the energy in the post-shock power law 
distribution of relativistic electrons that is radiated away. Here e# = (I? ps ) 2 /87re ps which for 
a sub-equipartition field is < 1, but for a field advected from the source and compressed by 
the shock it is (see above) ~ 64er(A7 /7)~ 2 which is ^> 1 unless the magnetization of the 
outflow is extremely low, a < 10 _2 (A7/7) 2 . Thus Y ~ e e e rad 10~ 2 o- _1 (A7/7) 3 < 1. For 
a shock-generated field that is a constant fraction of equipartition Y would also be small 
unless e rad e e /e B > 7/A7. 

The cooling break frequency thus scales as 

eB , si 2 [ (1 + Y)-h B 3/2 (^) 4 L^VAtej (equipartition) , 
7 o _^ oc ^ (All) 
K)- 2 (7- 3 / 2 (^V^r^VAtej (advected field) . 




Finally, we derive the fraction, e, of the total energy that is converted into internal 
energy during the collision between the two shells (in the limit a <C 1 and A7/7 <C 1). 
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Conservation of energy and momentum read 



7imi + 7 2 m 2 
7i/3imi + 7 2 /3 2 m 2 



7/ M , 
7//5/M , 



(A12) 
(A13) 



where 7/ = (1 — 0f)~ x ^ 2 is the final Lorentz factor, and M = m 1 + m 2 + E' mt /c 2 where i?( nt 
is the internal energy that was produced in the collision, as measured in the rest frame of 
the merged shell, while its value in the lab frame is E int = 7/-E; nt . One obtains (3f from the 
ratio of the two equations, and for 7 ^ 1 we have 



and therefore 
Thus we have 



l-(3 f 



7i(l - (3i)mi + 7 2 (1 - (3 2 )m 2 ^ mi/71 + m 2 Ay 2 



7imi + 7 2 m 2 



2(7imi + 7 2 m 2 ) 



If 



7i™i + 7 2 m 2 



e = 



1 - 



^1/71 + ^2/72 
7/(mi + m 2 ) 



(71ml + 7 2 m 2 )c 2 71ml + 7 2 m 2 

and since for A7/7 <C 1 (in addition to 7 3> 1), 



1 - 



1 + 



m x m 2 



[mi + m 2 y 



72 

7i 



7i 

72 



2(72i - 1) 



72 
7i 



7i 
72 



then 



m\m 2 



( A7 

2 (mi + m 2 ) 2 \ 7 



x 





2(i+x)=» v 7 ; 



(A14) 

(A15) 

-1/2 
(A16) 

(A17) 
(A18) 



where x = m 2 /mi (or alternatively mi/m 2 ) is the rest-mass ratio of the two shells. 
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Table 1. Data displayed in the top three panels of Figure [5] 



Time Interval 


Fluence (0.3-10 keV) 


Fluence (15-150 keV) 


Epeak (keV) 


Power Law index 


E iso (10 52 ergs) 


-13.4 - 18.0 




, 05+0.45 

0.ttO_ Q gg 


> 77.6 


1.61 ± 0.13 


> 3.74 


70.21 - 86.2 
86.2 - 102.88 
107.0 - 121.04 


1.35± 0.16 


„ 44 +0.35 
1 qn +1 ' 68 
0.66 ± 0.15 


55.5 ± 11.5 
45.0 ± 9.6 
9.8 ± 1.5 


1.00 ± 0.51 
1.38 ± 0.36 
1.50 ± 0.16 


1 36+ ' 55 
± - ou -0 .29 

1 76+ ' 74 

L - ' D _0.39 

0.64 ± 0.07 


121.04 - 159.21 
159.21 - 199.21 


58+ ' 06 
u -°°-0.07 

32+ 34 
u -° -0.09 


06+ 06 
< 0.0504 


3.9 ± 0.77 
0.5 ± 0.50 


1.51 ± 0.24 
1.90 ± 0.26 


0.58 ± 0.10 
0.36 ±0.07 



Note. — Fluence is in units 10 8 erg cm 2 



Table 2. Temporal decay slopes of flares 



Flare 


Peak 


Start t\ 


Rise time At (s) 


Decay index ct\ 


to for q b (s) 


Decay index aa 


prompt 




-13.4 


15.5 


5.4 ± 1.1 


-13.17 


2.34 ± 0.45 


prompt decay 


18.0 




2.1 ± 0.8 


-13.17 


1.16 ± 0.47 


1 


1 


70.2 


0.19 


16.3 ± 5.3 


70.21 


0.29 ± 0.11 


1 


2 


73.6 


2.01 


24.0 ± 4.1 


70.21 


2.14 ± 0.36 


1 


3 


78.8 


0.39 


23.9 ± 2.8 


70.21 


3.35 ± 0.39 


2 


1 


86.2 


1.77 


19.2 ± 2.5 


86.2 


1.33 ± 0.18 


2 


2 


95.5 


3.54 


23.7 ± 4.8 


86.2 


3.50 ± 0.71 


3 


1 


102.9 


2.54 


46.9 ± 17.7 


108.5 


3.77 ± 1.42 


3 


2 


118.3 


3.14 


41.2 ± 18.2 


108.5 


4.56 ± 2.01 


4 


1 


122.5 


11.2 


10.3 ± 5.1 


122.54 


0.96 ± 0.46 


4 


2 


136.2 


2.52 


10.7 ± 0.5 


122.54 


1.78 ± 0.08 


5 




160.7 


18.0 


17.7 ± 0.7 


160.71 


2.73 ± 0.12 


final prompt decay 


195.6 




6.34 ± 0.39 


160.71 


2.14 ± 0.13 


shallow afterglow 


323.8 




0.31 ± 0.17 


160.71 


0.24 ± 0.03 


steep afterglow 


3200 




1.24 ± 0.05 


160.71 


1.21 ± 0.03 



Note. — Flares are numbered as in the text and Figures [T] and [2] Times are with reference to the trigger 
time To and the definition of the rise time is given in the text. The decay index is derived through a fit to: 
R oc (t — to)~ a , where R is the photon event rate, t is the time, and to is defined as the trigger time To when 
deriving a\ (column 5), and as the start of the particular peak or flare when deriving «b (column 7). The 
values of to used in the derivation of c«b are shown in column 6. 
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-50 50 100 150 200 250 300 

T-TO (s) 



Fig. 1. — The BAT light curve is shown in four energy bands and the sum of all energy 
bands. The vertical bars indicate the start and end of the spacecraft slew to the burst 
location. Note that the count rate statistical errors are much larger before the slew than 
after the slew. This is because the burst was detected near the edge of the BAT field of view 
where only 27% of the detectors were illuminated. None of the apparent sharp structure in 
the prompt emission is statistically significant. The arrows in the bottom plot indicate the 
peak and the numbering of each of the flares. 
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Fig. 2. — The five flares are shown for both BAT (open squares) and XRT (crosses). The first 
XRT flare at T ~ 115 s is clearly detected in the BAT and there appears to BAT emission at 
the peak of the flare at T ~ 175 s. The arrows indicate the peak and the numbering of each 
of the flares. The time for which the spectral fits to the five flares were made are indicated 
by the bars at the bottom of the plot. 
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Fig. 3. — The XRT light curve showing the early time flares and three episodes of smoothly 
decaying afterglow emission. 
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Fig. 4. — The optical measurements (various symbols; right hand scale) and X-ray count 
rate (crosses; left hand scale) for the late time observations. The fits to the data discussed 
in the text are shown as solid lines. The optical symbols are defined as follows: V: filled 
triangles, White: open triangle s, B: open ci rcles, R: filled circle s, J: open square, I : fille d 
square. Credits: I and J band: ICobb (120061) : unlabled R band : lAsfandvarof et al. (2006); 



J2OO6 re- 



l abeled poi nts: a. Asfandyarof et al. ( 2006h. b. Pavlenko et al. ( 2006a ) 



d. Pavlenko et al 



(l2006bh .e. iRumyantsev et al.l 



c. 



Jakobsson et al. 



Jakobsson et al. 



d2006ch . 



Jakobsson et al.l (l2006dl ). All other points are UVOT measurements. U and Wl band 



upper limits are omitted for clarity. All upper limits are 2a. 
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Fig. 5. — The evolution of spectral fit and energetic properties of the flares. The fits are to: circles 
- BAT alone, triangles - BAT and XRT jointly. We were unable to fit a cut-off power law to the 
prompt emission so the points for the prompt emission show lower limits on Speak and E[ so . The top 
panel shows how -Speak (shown in the observer frame) changes from the prompt emission through 
the five flares. The second panel shows the evolution of the power-law index (a in the cut-off power 
law fit) across the flares. Both of the two top plots clearly show a hard to soft spectral evolution as 
the flares progress. The third panel shows the time evolution of the isotropic radiated energy over 
the 1-10 4 keV energy range, indicating the the flares become progressively less energetic. The data 
plotted in the top three panels is also given in Table [TJ The bottom panel shows hardness ratios for 
the burst. For the prompt emission and the first four flares the flux ratio S(50-100 keV)/S(25-50 
keV) is shown as circles. (The hardness ratio for the last flare is consistent with zero.) For the last 
three flares, we also show the flux ratio S(1.5-10 keV)/S(0.3-1.5 keV) (triangles). 
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Fig. 6. — The peak energy is plotted ag ainst the isotropic e nergy over the 1-10 4 keV energy 
range (the same ener gy range used b y lAmati et al.ll2002l ). The dashed line is the fit to 
E pmk -E im derived by iGhirlanda et al.l (12004 ) and the dot-dashed line is the fit derived by 



Amati et al.l (120021 ). The time ordering of the flares goes monotonically from highest E pea k 



to lowest. The prompt emission and the first two flares detected fall on the Amati relation 
while the last three flares fall below 
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Fig. 7. — The prompt emission (top), the BAT-detected flares (middle) and the XRT- 
detected flares (bottom) are shown on the same time scale. The solid lines on the plot show 
the best power law fit to each individual segment of the light curve and the changes in slope 
indicate the start, apex and end of each subpeak. The method for defining the intervals is 
given in the text and decay constants are listed in Table O One can note the sharp temporal 
features (subpeaks) of the flares. It is also clear that the X-ray flares are longer than the 
BAT flares, a feature also seen in lRomano et al.l (120061 ). and consistent with previous results 
showing that the duration of pulses in prompt emission are longer at low energy than at 
high energy. 
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Fig. 8.— Plot of At/t vs. AF/F. The flux ratio AF/F is derived as AF = (F(t 2 ) - F(ti)) 
and F = F(tx), where F(ti) and F(t 2 ) are the flux at the start and top of each peak, 
respectively. In this notation AF/ F = 1 means a doubling of the flux. This figure shows that 
for all peaks AF/F ~ 1 while At/t <C 1. The area to the right of th e dashed li n es ind icates 
the kinematically allowed region for afterglow variability derived by lloka et al. dgQoA The 



vertic al line indicates that refreshed shocks cannot make a bump with At < t/4. lloka et al. 
( 120051 ) also argue that ambient density fluctuations cannot make a bump in afterglow light 
curves larger than the limit indicated by the diagonal line. 
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Fig. 9. — This plot shows At (rise time to the peak) and At/t vs. time since the GRB trigger 
for each peak in each of the flares. Note that the plot is log-log. The quantities At and t 
are defined in the text. Although there is some scatter there is a general trend for both At 
and At/t to increase with t for the flares. A fit to At vs. time (top plot) gives a slope of 
4.0 ± 0.9. It is also quite clear that for all peaks in the flares, that At/t <C 1 



